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We study the r-modes and rotational "hybrid" modes of relativistic stars. As in Newtonian 
gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation 
of state. If the star and its perturbations obey the same one-parameter equation of state (as with 
isentropic stars), there exist no pure r-modes at all - no modes whose limit, for a star with zero 
angular velocity, is an axial-parity perturbation. Rotating stars of this kind similarly have no pure 
g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed 
pressure and density. 

In spherical stars of this kind, the r-modes and g-modes form a degenerate zero-frequency sub- 
space. We find that rotation splits the degeneracy to zeroth order in the star's angular velocity 
Q, and the resulting modes are generically hybrids, whose limit as f2 — > is a stationary current 
with both axial and polar parts. Because each mode has definite parity, its axial and polar parts 
have alternating values of I. We show that each mode belongs to one of two classes, axial-led or 
polar-led, depending on whether the spherical harmonic with the lowest value of I that contributes 
to its velocity field is axial or polar. Newtonian isentropic stars retain a vestigial set of purely axial 
modes (those with I = m); however, for relativistic isentropic stars we show that these modes must 
also be replaced by axial-led hybrids. We compute the post-Newtonian corrections to the I — m 
modes for uniform density stars. 

On the other hand, if the star is non-isentropic (or, more broadly, if the perturbed star obeys 
an equation of state that differs from that of the unperturbed star) the r-modes alone span the 
degenerate zero-frequency subspace of the spherical star. In Newtonian stars, this degeneracy is 
split only by the order Q,^ rotational corrections. However, when relativistic effects are included the 
degeneracy is again broken at zeroth order. We compute the r-modes of a non-isentropic, uniform 
density model to first post-Ne-wtonian order. 



I. INTRODUCTION 



The discovery that the r-modes in rotating stars are genericaUy unstable due to the emission of gravitational -waves 
has attracted a large amount of attention in the last t-wo years. The current models suggest that the r-mode 
instability may cause a ne-wly born neutron star to spin do-wn to a fraction of the Kepler frequency (that provides the 
limit of dynamical stability) in the first fe-w months of its existence . Since a considerable amount of gravitational 
radiation is generated in the process, the r-modes provide a promising source for the generation of gravitational -wave 
interferometers that are currently under construction ||^ . It is also speculated that the instability associated -with the 
r-modes may be relevant for older neutron stars in accreting systems 1^,0. 



*Email address: lockitch@gravity.phys.psu.edu 
^Email address: na@maths.soton.ac.uk 
■'■Email address: friedman@uwm.edu 



1 



Since the instability was first discovered, and its potential astrophysical relevance was appreciated, there have been 
many attempts to improve on the detailed physics incorporated in the models. This effort leads to difficult questions 
re gard ing, for example, neutron star superfluidity the interplay between the magnetic field and fluid pulsations 
|p|-|rT|, and the formation of a solid crust as a young neutron star cools [p^-p^. These, and several other, issues must 
be addressed before the true astrophysical relevance of the r- modes can be assessed. Our understanding of the r-mode 
instability, however, is based almost entirely on Newtonian calculations, and it is important to compute these modes 
in a relativistic context, where instability growth times may differ significantly from Newtonian-based estimates. (The 
closely related instability of f-modes of rapidly rotating stars is sharply strengthened by relativistic effects; see |l9| ] 
for a review.) ^ 

The purpose of the present investigation is to understand how general relativity affects the properties of the r-modes. 
In order to address this issue we first need to discuss the general nature of modes of rotating stars. 

The spherical symmetry of a non-rotating star implies that its perturbations can be divided into two classes, polar or 
axial, according to their behaviour under parity. Where polar tensor fields on a 2-sphere can be constructed from the 
scalars y";™ and their gradients VYJ™ (and the metric on a 2-sphere), axial fields involve the pseudo- vector f x VYJ™, 
and their behavior under parity is opposite to that of Ij™. That is, axial perturbations of odd I are invariant under 
parity, and axial perturbations with even I change sign. Because a rotating star is also invariant under parity, its 
perturbations may also be divided into distinct parity eigenstates. If a mode varies continuously along a sequence 
of equilibrium configurations that starts with a spherical star and continues along a path of increasing rotation, the 
mode will be called axial if it is axial for the spherical star. Its parity cannot change along the sequence, but I is 
well-defined only for modes of the spherical configuration. 

It is useful to subdivide stellar pulsation modes according to the physics dominating their behaviour. This classifi- 
cation was first developed by Cowling pl| for the polar perturbations of Newtonian polytropic models. The p-modes 
of spherical models are polar-parity modes having pressure as their dominant restoring force|^ They typically have 
large pressure and density perturbations and high frequencies (higher than a few kilohertz for neutron stars). The g- 
modes are polar-parity modes that are chiefly restored by gravity. They typically have very small pressure and density 
perturbations and low frequencies. Indeed, for spherical isentropic stars, which are marginally stable to convection, 
the g-modes are all zero-frequency and have vanishing perturbed pressure and density. Similarly, all axial-parity 
perturbations of non-rotating perfect fluid models have zero frequency. The perturbed pressure and density as well 
as the radial component of the fluid velocity all vanish for axial perturbations; being rotational scalars, they must 
have polar parity. Thus, the axial perturbations of a spherical star are simply stationary horizontal fluid currents. 
This Newtonian picture of stellar pulsation is readily generalised to the relativistic case. The only difference is that 
the various modes will now generate gravitational waves. This means that they are no longer "normal modes" , but 
satisfy outgoing wave boundary conditions at spatial infinity. Furthermore, one can identify an additional class of 
such outgoing modes in relativistic stars. Like the modes of black holes, these modes are essentially associated with 
the dynamical spacetime geometry and have been termed w-modes, or gravitational- wave modes p^ . For a general 
discussion of the oscillations of relativistic stars we refer the reader to the recent review article by Kokkotas and 
Schmidt ||. 

In general, the classification of modes is relevant also for rotating stars, even though the character of the various 
modes may be significantly affected by rotation. In particular, rotation imparts a finite frequency to the zero- 
frequency perturbations of spherical stars. Because these modes are restored by the Coriolis force, their frequencies 
are proportional to the star's angular velocity, f2. In fluid mechanics, such modes are generally known as inertial 
modes [^,^^. In non-isentropic stars these rotationally restored modes all have axial parity (the polar g-modes 
are nondegenerate already in a spherical non-isentropic star because of internal entropy gradients). In astrophysics, 
these modes were first studied in Newtonian gravity by Papaloizou and Pringle |26| , who called them r-modes because 
of their similarity to the Rossby waves of terrestrial meteorology. In isentropic stars, however, the space of zero 
frequency modes of the spherical model includes the polar g-modes in addition to the axial r-modes. This large 
degenerate subspace of zero-frequency modes is split by rotation to zeroth order in the angular velocity, and the 
rotationally restored (inertial) modes of isentropic stars are generically hybrids whose spherical limits are mixtures of 
axial and polar perturbations. This has been shown in Newtonian gravity by Lockitch and Friedman p7| (see also 
1^,^). In order to distinguish between the two classes of inertial modes we refer to modes which become purely 
axial in the spherical limit as r-modes, while modes that limit to a mixed parity state are called rotational hybrid 



^There are as yet no fully relativistic calculations of other pulsation modes (like the f-mode) of rapidly rotating stars, except 
in the Cowling approximation ]2c| ]. 
■^The lowest p-mode for each value of / and m is termed an f-mode or fundamental mode. 
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modes. This is a natural nomenclature given the standard distinction between axial and polar modes in relativistic 
studies of spherical stars. 

Attempts to study the r- modes of rotating relativistic stars were not made until rather recently [|l 30-3^. In fact, 



the present investigation is the first study of this problem that puts all its different facets in the proper context. In 
particular, we prove that (apart from a set of stationary dipole modes) rotating relativistic isentropic stars have no 
pure r- modes (modes whose limit for a spherical star is purely axial) . This is in contrast with the isentropic Newtonian 
stars which retain a vestigial set of purely axial modes (those having spherical harmonic indices I — to). Instead, the 
Newtonian r-modes with I = m > 2 acquire relativistic corrections with both axial and polar parity to become discrete 
hybrid modes of the corresponding relativistic models. We compute these corrections for slowly rotating isentropic 
stars to first post-Newtonian order. 

For non- isentropic relativistic stars the situation is somewhat different. In the slow- motion approximation in which 
they have so far been studied, non-isentropic stars have, remarkably, a continuous spectrum. Kojima ]30| ] has shown 
that purely axial modes would be described by a single, second-order ODE for the modes' radial behavior. He then 
argues that the continuous spectrum is implied by the fact that eigenvalue problem is singular (the coefficient of 
the highest derivative term of the equation vanishes at some value of the radial coordinate). This claim has been 
made mathematically precise by Beyer and Kokkotas As the latter authors point out, the continuous spectrum 
may be an artifact of the vanishing of the imaginary part of the frequency in the slow rotation limit. (Or, more 
broadly, it may be an artifact of the slow rotation approximation itself.) In this paper we show that, in addition to 
the continuous spectrum, certain discrete modes also exist as solutions to Kojima's equation. These modes are the 
relativistic analogue to the Newtonian r- modes in non-isentropic stars. We compute these modes for slowly rotating 
non-isentropic stars. 

In a complementary study of the relativistic r-modes, Kojima and Hosonuma [ p4| have recently derived the order 
rotational corrections to Kojima's equation. Working in the time domain, they derive a set of evolution equations for 
an axial perturbation and its lowest order polar and axial corrections. Direct numerical evolution of these equations 
(with appropriate initial data) would provide a useful comparison with our results on the modes of non-isentropic 
relativistic stars. 

When does one need to take into account the departure of a neutron star from isentropy in computing rotational 
modes? Because of bulk viscosity, a gravitational- wave driven instability is unlikely to set in above about 10^'-' K. This 
is well below the Fermi temperature of the star's baryons, and the departure from isentropy appears to be dominated 
by composition gradients in the crust and interior. These have been discussed in the context of g-modes of spherical 
stars by Finn Q and by Reisenegger and Goldreich |3^j3^ and for rotating stars by Lai Q . Because the timescale 
of perturbations is too slow to allow the beta- and inverse beta-decays that would allow a displaced fluid element to 
adjust its composition to that of the surrounding star, Ap/Ap is greater than d\ogp/d\og p by a factor 1 -I- ^x, where 
X — n-p/n « 6 X 10~^p/pnuciear IS the local ratio of protons to baryons. This leads in the star's interior to g-mode 
frequencies limited by the Brunt- Vaisala-frequency, 

q \ 1/2 

^ .^nn.-V P ^1/2. 
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when a crust is present, crustal g-modes have comparable frequencies. The g-mode frequencies of spherical stars are 
then of order 100-200 Hz; when this is smaller than the frequencies of the rotationally restored modes of the isentropic 
models, one expects the isentropic approximation to be valid. 

The plan of the paper is as follows. We begin, in Section |l[ with a brief review of the Eulerian and Lagrangian 



perturbation formalisms, both of which are used in the paper. In Section III, we consider the time-independent 
perturbations of spherical relativistic stars and prove that the subspace of nonradial zero-frequency modes is spanned 
by the r- and g-modes in isentropic models, but by the r-modes alone in non-isentropic models. Because of this 
difference, the character of th e mo de spectrum in rotating isentropic models differs considerably from that of non- 



isentropic models. In Section IV A, we consider rotating non-isentropic stars and argue that the problem of finding 



their r-modes is well-defined. In Section IV B, we consider the isentropic case and derive a set of perturbation 
equations whose structure parallels the corresponding Newtonian equations of Lockitch and Friedman [p7|p|. This 
similarity between the Newtonian and relativistic equations leads to an identical structure of the mode spectrum and 
to a parallel theorem that every non-radial mode is either an axial-led or polar-led hybrid (the result has so far been 
proven only for slowly rotating relativistic stars). We consider the relativistic r-mode solutions of isentropic stars. 



^We will refer to equations from Ref. |2^] by the equation number with the prefix "LF". For example, Eq. (LF,25) will mean 
equation (25) from Ref. p.q. 
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finding that the zero-frequency dipole {1 = 1) solutions are the only purely axial solutions allowed. In other words, 
there are no nonstationary modes in isentropic relativistic stars whose limit as — > is a pure axial perturbation. In 
particular, the Newtonian r-modes having I = m > 2 do not exist in isentropic relativistic stars and must be replaced 
by axial- led hybrid modes . Thi s section concludes with a discussion of the boundary conditions appropriate to the 
relativistic modes (Section IV C). Finally, in Section^ we construct the post-Newtonian corrections to the well-known 
Newtonian r-modes in uniform density stars, both isentropic and non-isentropic. Some of the detailed equations, as 
well as the proof of the theorem regarding the isentropic mode spectrum, are presented in Appendices A-C. We use 
geometrized units {G = c = 1) throughout the paper. 



II. EULERIAN AND LAGRANGIAN PERTURBATIONS 



In general relativity, a complete description of a self-gravitating perfect fluid configuration is provided by a spacetime 
with metric ga/i, sourced by an energy-momentum tensor. 



Taf3 = (e + p)UaU0 + pgafl, 

where the fluid 4-velocity is a unit timelike vector field, 

u"Ury = —1, 



(2.1) 



(2.2) 



and e and p are, respectively, the total energy density and pressure of the fluid as measured by an observer moving 
with 4-velocity u". The metric and fluid variables satisfy an equation of state, 



e = s), 

with s the entropy per baryon, as well as the Einstein equation, 



(2.3) 



(2.4) 



An equilibrium stellar model is a stationary solution {gap, u" , e,p) to these equations. In this paper we will consider 
only equilibrium models obeying a barotropic (one-parameter) equation of state, 



(2.5) 



because this accurately models the equilibrium configuration of a neutron star. 

Adiabatic perturbations of such a star may be studied using either the Eulerian or the Lagrangian perturbation 
formalism 40 1. An Eulerian perturbation may be described in terms of a smooth family, [gc(/3(A), ?I"(A), e(A),p(A)], 
of solutions to the exact equations (2.2)-(2.4) that coincides with the equilibrium solution at A = 0, 

[5a/3(0),M"(0),e-(0),p(0)] = (g„0,u",e,p). 

Then the Eulerian change 5Q in a quantity Q may be defined (to linear order in A) as, 

dQ 



SQ 



(2.6) 



A=0 



Thus, an Eulerian perturbation is simply a change {hai3,Su°',Se,6p) in the equilibrium configuration at a particular 
point in spacetime (where we have written the change in the metric as haf} = 5 gap)- These must satisfy the perturbed 
Einstein equation 5G£ = SttSTJ^, together with an equa tion of state relating 5e and 5p that may, in general, differ 
from that of the equilibrium configuration (see Eq. ( 2.13 ) below). 

In the Lagrangian perturbation formalism | |39| , p0| , on the other hand, perturbed quantities are expressed in terms 
of the Eulerian change in the metric, hap and a Lagrangian displacement vector which connects fiuid elements in 
the equilibrium star to the corresponding elements in the perturbed star. The Lagrangian change AQ in a quantity 
Q is related to its Eulerian change SQ by 



AQ = 6Q + £^Q, 



(2.7) 



with the Lie derivative along ^° 
The identities, 
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A5„;3 = /i„/3 + 2V(„e/3) (2.8) 

A£„/3^5 = ^EapfS ff^^A^^^ (2.9) 

then allow one to express the fluid perturbation in terms of haf3 and 

Au" = iw^u^uTAg^^ (2.10) 

L ip e + p n 



where Fi is the adiabatic index, n is the baryon density and q"^ = g"^ + u^u^. Using Eqs. (2.7)-(2.11), it is 
straightforward to express the corresponding Eulerian changes also in terms of haf3 and e.g., 

Su°' = q°'p£ui,'^ + lu'^u^u'^hp^. (2.12) 



For an adiabatic perturbation of a barotropic equilibrium model, Eqs. (2/7) and ( ^.11 ) imply that the Eulerian 
changes in the pressure and energy density are related by 

^ = (2.13) 

Tip [e+p) 

where we have introduced the Schwarzschild discriminant, 

Ac. = V„e - ^Vop, (2.14) 

(e+p) Vip 

which governs convective stability in the star. In general, the adiabatic index Fi need not be the function 

F^ii±^^ (2.15) 
p at 

associated with the equilibrium equation of state. In terms of this function we have 

^o.= {^- ^V„p. (2.16) 

We will call a model isentropic if the perturbed configuration satisfies the same barotropic equation of state as the 
unperturbed configuration. In this case, Fi = F and the Schwarzschild discriminant vanishes identically. Such stars are 
marginally stable to convection. In this paper we study low-frequency pulsation modes of slowly rotating relativistic 
stars. We consider both isentropic and non-isentropic models. 

III. STATIONARY PERTURBATIONS OF SPHERICAL STARS 

The equilibrium of a spherical perfect fluid star is described by a static, spherically symmetric spacetime with 
metric gap of the form. 



and the energy- momentum tensor (2.1) with the fluid four- velocity given by 

u" = e'^t" . (3.2) 

Here t" = {dt)" is the timelike Killing vector of the spacetime. 

For barotropic stars, the pressure and energy density are related by an equation of state of form 

p = p{e). (3.3) 

In addition to this, the various quantities must satisfy the Einstein equation, Gafj = SnTap, which leads to the 
standard TOV equations. 
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dp _ (e + p) (M + Anr^p) 
dr " r{r - 2M) 

dM 



(3.4) 



Airr^e , (3.5) 
dr 



and 

dv 1 dp 



(3.6) 



dr (e +p) dr 
where 

Af(r) = ir(l - e-2^) . (3.7) 

Our main focus in this study is on the low-frequency osciUations, corresponding to r-modes and hybrid modes, 
of slowly rotating stars. As in Newtonian theory, we expect these modes to limit to stationary perturbations of a 
spherical star as the rotation rate goes to zero. In other words, we are interested in the space of zero-frequency modes; 
the linearized, time- independent perturbations of the static equilibrium. As in the Newtonian case |2^], we find that 
this zero-frequency subspace is spanned by two classes of perturbations. To identify these classes explicitly, we must 
examine the equations governing the perturbed configuration. 

Using the Eulerian formalism, we express the perturbed configuration in terms of the set (ft.Q/3, Su", Se, Sp), satisfying 
the perturbed Einstein equation SG^ = SttSTJ^ , together with an equation of state relating Se and 5p. 

The perturbed Einstein tensor is given by 

+ 2i?J/i/ + {V^V^Kp - V^VT/i - R''f'Kp)5£^ (3.8) 

where h = g"^haf3, Vq, is the covariant derivative associated with the equilibrium metric and 

r£ = 8tt{T/ - ^T6f) = 8TT [{e+p)uau^ + i(e -p)^/] (3.9) 

is the equilibrium Ricci tensor. The perturbed energy-momentum tensor is given by 

5T/ = {6e + 5p)u„uf^ + 5p5£ + (e + p)6u^uf^ + (e + p)uju'^. (3.10) 

Following Thorne and Campolattaro | |4l[ |, we expand our perturbed variables in scalar, vector and tensor spherical 
harmonics. The perturbed energy density and pressure are scalars and therefore must have polar parity 

5e - 6e{r)Yr, (3.11) 



5p = 5p{r)Yr. 

The perturbed 4-velocity for a polar-parity mode can be written 



(3.12) 



(3.13) 



while that of an axial-parity mode can be written 

<5< = -U{r)eS^~''h''P''^V pY{^u^V s r. (3.14) 

(We have chosen the exact form of these expressions for later convenience.) 

To simplify the form of the metric perturbation we will again follow Thorne and Campolattaro |41| and work in 
the Regge- Wheeler E2| gauge. The metric perturbation for a polar-parity mode can then be written 
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symm 






Hi{r) 










^K{r) 

r'^sin^eK{r) 



(3.15) 



while that of an axial-parity mode can be written 



/,o(r)(-i)9^r," 




symm symm 
symm symm 







hi{r) smOdeYl' 





(3.16) 



The Regge- Wheeler gauge is unique for perturbations having spherical harmonic index I > 2. However, when I — 1 
or I = 0, there remain additional gauge degrees of freedom^ In addition, the components of the perturbed Einstein 
equation acquire a slightly different form in each of the three cases. (Campolattaro and Thorne discuss the 
difference between the I > 2 and I = 1 cases.) 

We have d erive d the c omponents of t he perturbed Einstein equation using the Maple tensor packa ge^ b y substituting 



and ( |3.10D (making liberal use of the equilibrium equations (3.4) through ( |3.7| ) 
The resulting set of equations for the case I > 2 are equivalent to those presented by 



expressions (3.11)-(3.16) into Eqs. (3. 
to simplify the expressions) 

Thorne and Campolattaro |41| ] upon specializing their equations to the case of stationary perturbations and making 
the necessary changes of notationrl Similarly, the set of equations for the case I — 1 are equivalent to those presented 
by Campolattaro and Thorne [ [43[ . For completeness, the equations governing stationary perturbations of spherical 
stars are given in Appendix A. 



A. Decomposition of the zero-frequency subspace. 



By inspection of the three sets of perturbation equations given in Appendix A, it is evident that they decouple into 
two independent classes. We find that any solution 



{Ho,Hi,H2,K,ho,W,V,U,5e,6p), 



(3.17) 



to the equations governing the time-independent perturbations of a static, spherical star is a superposition of (i) a 
solution 



iO,Hi,0,0,ho,W,V,U,0,0) 



to Eqs. (|A6|)-(|A§) or ( \A2l\ ) and (ii) a solution 

{Ho,0,H2,K,0,0,0,0,Se,5p) 



(3.18) 



(3.19) 



to Eqs. (Q-(^), (|aTT1)-(|aT^) or (1aT8|)-([A20|). 

For solutions of type (ii), the vanishing of the (tr), {t9) and (tip) components of the perturbed metric in our 
coordinate system implies that these solutions are static. If one assumes the linearization stability]] of these solutions, 
i.e., that any solution to the static perturbation equations is tangent to a family of exact static solutions, then the 



^Letting cab be the metric on a two-sphere with eab and Da the associated volume element and covariant derivative, 
respectively, one finds the following: When Z > 2 the polar tensors DaDbYi™ and CABYf" are linearly independent, but when 
1 = 1, they coincide. In addition, the axial tensor t^^^ Dc)DbY{^ vanishes identically for 1 = 1 and, of course, DaY{^ vanishes 
for I = 0. 

^ www. maplesoft . com 

^In particular, their fluid variables (denoted by the subscript tc) are related to ours as follows: UTciryt) = Ut, WTc{r,t) = 
—re^Wt, VTcif^t) = Vt and St/{e -\- p) = 5p/Vip = —{K -\- |i/o). Their equilibrium metric has the opposite signature and 
differs in the definitions of the metric potentials vtc = \i' and Arc = 

''We are aware of a proof of this linearization stability property under assumptions on the equation of state that are satisfied 
by uniform density stars, but would not allow polytropes mM. 



7 



theorem that any static self-gravitating perfect fluid is spherical implies that any solution of type (ii) is simply a 
neighboring spherical equilibrium. 

Thus, under the assumption of linearization stability we have shown that all stationary non-radial {I > 0) pertur- 
bations of a spherical star have 

Ho = H2 = K = Se = 5p = 



and satisfy Eqs. (|A6|)-(|A8|). That is, 



= H, + l?Si±i^e^Vw^, 



1(1 + 1) 



= e-^"--^) 



+ 167r(e + p)e^^V, 



(3.20) 
(3.21) 



(i/' + X')h'o + 



ho - -{u' + X')U 



(3.22) 



where a prime denotes a derivative with respect to r. Observe that if we use Eq. ( |3.20 ) to eliminate Hi{r) from Eq. 
(3.21) we obtain 



V 



l{l + l){e + p) 



[(e-fp)e'^+^rM^]'. 



(3.23) 



This equation is clearly the generalization to relativistic stars of the conservation of mass equation in Newtonian 
gravity, Eq. (LF,13). The other two equations relate the two dynamical degrees of freedom of the spacetime metric 
to the perturbation of the fluid 4- velocity and vanish in the Newtonian limit. 

These perturbations must be regular everywhere and satisfy the boundary condition that the Lagrangian change in 
the pressure vanishes at the surface of the star, r — R. We show in Section IV C below that this boundary condition 
requires only 



W{R) = 0. 



(3.24) 



leaving W{r) and U{r) otherwise undetermined. If W{r) and U{r) are specified, then the functions Hi{r), ho{r) and 
V{r) are determined by the above equations. The solutions for the metric variables are subjec t to m atching conditions 
to the solutions in the exterior spacetime, which must also be regular at infinity, see Section IV C . 

Finally, we consid er the equation of state of the perturbed star. We have seen that for an adiabatic oscillation of a 
barotropic star Eq. (2.11) implies that the perturbed pressure and energy density are related by 



(5e 



5p 

Tip (e+p) 



(e + p) Tip 



for some adiabatic index ri(r). 

The Lagrangian displacement vector ^" is related to our perturbation variables by Eq. ( ^.12 ), 



Thus, we have 



or [taking the initial displacement (at f = 0) to be zero] 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



For the class of perturbations under consideration, we have seen that 8p ~ 8e ~ Thus Eqs. ( 3.25 ) and ( 3.2g 
require that 



_(e+p) Tip_ 



0. 



(3.29) 
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For axial-parity perturbations this equation is automatically satisfied, since Suj^ has no r-component, cf. Eq. (3.14) 
In other words, a spherical barotropic star always admits a set of a xial z ero-frequency modes (the r- modes). 
For polar-parity perturbations, 5up = e~''W{r)/r ^ 0, and Eq. ( 3.29 ) will be satisfied if and only if 



r,(r).r(.)=(i±^^. (3.30) 

p ae 

Thus, we see that a spherical barotropic star admits a class of zero-frequency modes (the g-modes) if and only if the 
perturbed star obeys the same one-parameter equation of state as the equilibrium star. Again, we call such stars 
isentropic, because isentropic models and their adiabatic perturbations obey the same one-parameter equation of 
state. 

That all axial-parity fluid perturbations of a spherical relativistic star are time-independent was flrst shown by 
Thorne and Campolattaro [}4l| . The time- independent g-modes in spherical, isentropic, relativistic stars were found 
by Thorne |4|. 

Summarizing our results, we have shown that a spherical barotropic star always admits a class of zero- frequency r- 
modes (stationary fluid currents with axial parity) ; but admits zero- frequency g-modes (stationary fluid currents with 
polar parity) if and only if the star is isentropic. Conversely, the zero-frequency subspace of non-radial perturbations 
of a non-isentropic spherical star is spanned by the r-modes only; while the zero-frequency subspace of non-radial 
perturbations of a spherical isentropic star is spanned by the r- and g-modes - that is, by convective fluid motions 
having both axial and polar parity and with vanishing perturbed pressure and density. Being stationary, these r- and 
g-modes do not couple to gravitational radiation, although the r-modes do induce a nontrivial metric perturbation 
{hte,htip 7^ 0) in the spacetime exterior to the star (frame-dragging). One would expect this large subspace of modes, 
which is degenerate at zero-frequency, to be split by rotation, as it is in Newtonian stars. Our aim is to investigate 
this issue in more detail, and we will begin by considering perturbations of slowly rotating relativistic stars. 



IV. PERTURBATIONS OF SLOWLY ROTATING STARS 



The equilibrium of a perfect fluid star that is rotating slowly with uniform angular velocity Q is described ^6 47 
by a stationary, axisymmetric spacetime with metric, gap, of the form 

ds^ = -e^^'^^'Ut^ + e^^^^'Ur^ + r^dO^ + r^sin^Odip^ - 2uj{r)r^ sin^Odtdip (4.1) 



(accurate to order Jl). The energy-momentum tensor follows from (2.1) and the fluid 4-velocity to order Vl 



u° = e-^ie +nip°') (4.2) 

Here t" = {dt)" and ip" — (dip)" are, respectively, the timelike and rotational Killing vectors of the spacetime. 

That the star is rotating slowly corresponds to the assumption that Q is small compared to the Kepler velocity, 
^Ik ~ \/ M/R^, the angular velocity at which the star is dynamically unstable to mass shedding at its equator. In 
particular, we neglect all quantities of order il^ or higher. To order fl the star retains its spheri cal sh ape , because the 



centrifugal deformation of its figure is an order fl^ effect This means that the equations (3^)- (3/7) governing a 



spherical star remain relevant also for a slowly rotating equilibrium configuration. In addition we need to solve an 



equation |46| that determines the new metric function uj{r) in terms of the spherical metric functions iy{r) and A(r), 



^^"^^^ W(-^)?V-f$^ + ^)^-0 (4.3) 



r** dr \ dr J r \ dr dr 

where 

Cj{r) = VL^Lo{r). (4.4) 

This new metric variable is a quantity of order SI that governs the dragging of inertial frames induced by the rotation 
of the star [Q. Apart from the frame-dragging effect, however, the spacetime is unchanged from the spherical 
configuration. Outside the star, Eq. ( |4.3| ) has the solution. 



where J is the angular momentum of the spacetime. This relation can be used to provide boundar y co nditions for ui 
(and its derivative) at the surface of the star in terms of and J. Specifically, the solution to Eq. (O) is normalized 
by requiring that 



u{R) + -Ru>'{R) 



(4.6) 



where R is the radius of the star. 

Note that ojc and u) satisfy the inequalities < ujc < OJ < ^ (where an index c denotes the value at the centre of 
the star). This means that < w < 57 — ojc, and that fl, uj and w are positive for all values of r. Defining a rescaled 
variable a) = (D/f2, we have Cjc = uic/^ < < 1. Then, to linear order in f2, a single integration of (4.3) suffices to 
determine the frame dragging for all 17 and a specific stellar model (a given equation of state and, say, the central 
density) . 

We now consider the nonradial perturbations of these slowly rotating equilibrium models to linear order in Q. Since 
the equilibrium spacetime is stationary and axisymmetric, we may decompose our perturbations into modes of the 
formUe*^'^*^'"'''^ The perturbation equations have been written down in the Eulerian formalism by Kojima ji^, but 
we will find it convenient to work also in the Lagrangian formalism. We therefore begin by expanding the perturbed 
density and pressure, the displacement vector and the metric perturbations hap in tensor spherical harmonics. 

The Eulerian change in the density and pressure may be written as 



56=^ <5e,(r)l^^ 



(4.7) 



and 



(4.8) 



respectively. 

The Lagrangian displacement vector can be written 



^ EE — V<^ -VF,(r)Y;™r" + V,(r)V"Y,'" - iUi{r)P'^^£t''^'^'^VpY/''V^tVsr\t 

l=m ^ ' 



(4.9) 



where we have defined, 



and introduced the comoving frequency, 



(4.10) 



(4.11) 



The exact form of expression (^^) has been chosen for later convenien ce. In particular, we have chosen a gauge in 
which = 0. Note also the chosen relative phase between the terms in (4.9) with polar parity (those with coefficients 
Wi and Vi) and the terms with axial parity (those with coefficients [//). 

Working in the Regge- Wheeler gauge, we express our metric perturbation as 



Hu{r)yi 



2\ 



syinm 
syinm 



H2Ar)e 

symm 
syinm 



Yr h^Ar)i±g)Yr ihuir) ^v^edeY- 
r^Ki{r)Yi"' 

r^sin^eKi{r)Y/'' 



(4.12) 



Again, note the choice of phase between the polar-parity components (those with coefficients Hq^i, Hi^i, H2,i and Ki) 
and the axial-parity components (those with coefficients /iq.; and /ii,;)- 



We will always choose m > since the complex conjugate of an m < mode with real frequency cr is an m > mode with 
frequency —a. Note that a is the frequency measured by an inertial observer at infinity. 
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The use of the Lagrangian formaUsm introduces additional gauge freedom into the problem of stellar perturbations. 
This freedom is associated with a class of trivial displacements that leave all physical quantities invariant p9| , ^ . 
One eliminates this gauge freedom by restricting attention to the "canonical" displacements — those that conserve 
vorticity in constant entropy surfaces [|T]j3^. This conservation law, known as Ertel's theorem, is essentially the curl 
of the perturbed Euler equation and in general relativity has the form , 

A j£„co„0-^V[„nV^]pl =0 (4.13) 



or. 



where 



2 r 

£u^i^o.p = ^ {V[<,AnV/3]p + V[„nV^]Ap} , (4.14) 



u:^0 = 2V[„ (4.15) 



i s the relativistic vorticity. For our slowly rotating equilibrium star, Eq. (4.14) can be written using Eqs. (2.11) and 
(^T6|) as. 



zKfle-^Aw^^ = -v4,V[„rV^]Ap, (4.16) 



since Aa = ArV aX, cf. ( 2.161 ). Note that the three spatial components of the perturbed vorticity are not independent, 



being related by the identity 

V[„Ac^0^]=O. (4.17) 

We seek those modes that in the limit f2 — > belong to the zero-frequency subspace considered in the previous 
section. We have shown that such modes must have axial parity in non-isentropic stars, but may be either polar or 
axial in the iscntropic case. We will, therefore, require that our perturbation variables obey an ordering in powers of 

that reflects this spherical limit, 

Ui,h^,i ^ 0(1) 



1,1 



0(1) 

0(^(2) 



isentropic stars 
non-isentropic stars 



HQ,i,H2,i,Ki,hi,i,Sei,Spi,cT - 0{fl) 



(4.18) 



In addition to the new equations that arise at order f2, the zeroth order quantities must obey the zeroth order 
perturbation equations ( [3.2C ), ( 3.22| ) and ( 3.23| ), for all I. The degeneracy of the zero- frequency modes will be split 
at zeroth order if there is a subset of the 0(17) equations that involves only the 0(1) variables. While this occurs in 
Newtonian gravity only for isentropic stars p7|, in general relativity it occurs also for non-isentropic stars. 



A. The non-isentropic case 



In a search for the relativistic r-modes, Kojima has recently applied his general perturbation equations p^ ] 
to the case of a mode whose spherical limit is purely axial. Accordingly, he assumes an ordering of his perturbation 
variables in powers of il that agrees with our non-isentropic ordering (although he does not distin guish between the 
isentropic and non-isentropic cases). Kojima then finds that the zeroth order equation, Eq. (3.22), is joined at order 
f2 by an additional pair of equations, which can be written. 



l{l + l)<.i{(j + muj)e' 



-2v 



h'n 



O.l 



2imLLj'e ^^h^i = 



(4.19) 
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and 



im [iGirip + e) 



e "1,; 

l&TTmui{p + ey^e-'^^Ui = . (4.20) 
These two equations can be combined to give a second relation between the zeroth order variables /io,i and J7;, 



2M 



hi., 



im I iOTTip + e)r uje 

„2-2v 



2- —2^7 o / — 2i/— 2A 7 1 / — 21/— ^at/ I 

' ho^i — zroj e ho^i + ui r e "o ; I 



mil 



2mLU 



1(1 + 1) 



Ui 







(4.21) 



Kojima derived this equation from the perturbed Einstein equation, but as we will see in Section [VB the e quation can 
be writte n do wn immediately in the Lagrangian formalism as one of the spatial components of Eq. ( 4.16 ), Awgip = 
[cf. Eq. ( 4.53| )]. Inspection of the Eulerian equations (as, fo r exa mple, given by Kojima |4^) appears to suggest that 
there is no other equation in addition to Eqs. ( 0.22| ) and ( 4.21 ) that involves only the zeroth order axial variables 
/io,i and Ui. However, after a closer study we find that for iscntropic stars there is, in fact, a third such equation, 
implying that the system is overdetermined. While the existence of this third equation is obscured by the Eulerian 
formalism, it arises naturally in the Lagrangian framework as the other independent spatial component of Eq. (4.16). 



In non-isentropic stars, this equation couples the 0(1) variables occuring on the LHS to the 0{fl) variables (such as 
the perturbed pressure and density) appearing on the RHS. However, for isentropic stars the R HS vanis hes identically 
(since Ar = 0) and the third equation relating only the 0(1) variables emerges [cf. Eq. ( 4.54| ) or (4.55)]. 

Hence, for isentropic stars the assumption that the mode is purely axial as 17 ^ leads to an overdetermined 
system of equat ions. The appropriate spherical limit is therefore one that also includes the polar variables Wi, Vi 
and as in (4.18). For non-isentropic stars, on the other hand, the r-mode assumption appears to be consistent. 
Combining Eq. ( [4.2lD with Eq. ( |3.22| ) gives Kojima's "master" equation for /iq,;. 



mf2 



2mLU 










'1(1 + 1) AM 






e 


dr 







+ 167r(p + e)((T + mn)hQj = . 



(4.22) 



Kojima used this equation to argue that the r-mode spectrum of a relativistic star is continuous. The conclusion 
that the equation supports a continuous spectrum was shown with more mathematical rigour by Beyer and Kokkotas 
[ pl| . Basically, the continuous spectrum arises because (4.22) is a singular eigenvalue problem; the combination 
CT + mil — 2mCj/l{l + 1) may have a zero in the interval r e [0, oo]. It is interesting to ask whether the presence of a 
continuous part of the spectrum is physical or whether it is an artifact of the approximations we have introduced. That 
the latter may be the case can be argued for in the following way. To leading order in the slow-rotation expansion the 
mode frequency cr is a real valued quantity, but at higher orders it must h ave a lso an imaginary part (corresponding 
to dissipation due to gravitational wave emission). If we were to consider (4.22) for complex frequencies, the problem 
will be regular and there will likely be no continuous spectrum. The possible presence of a continuous spectrum is 
an interesting issue that should be investigated in more detail, but it is not the fo cus of the present study. What we 
want to emphasize here is that two important questions regarding equation (4.22) have not yet been answered. First 
of all, it has not been shown that the problem is well-defined. As we have already s tated , one can show that the 
system of equations is overdetermined for isentropic stars. This means that, equation (4.22) can only be relevant for 
non-isentropic stars. But in order to show that the equation is, indeed, relevant we must demonstrate that all other 
perturbation variables are uniquely specified given a solution for /iq,; from Eq. (4.22). Given the relative comp lexity 
of the perturbed Einstein equation, this is not a trivial task. Secondly, we need to investigate whether Eq. ( 4.22 ) 
allows distinct mode solutions in addition to its continuous spectrum. After all, the true relativistic analogue to a 
Newtonian r-mode ought to be a distinct mode with a well-defined eigenfunction. 

We address the first issue by considering the perturbation equations that ari se in th e Eu lerian formalism, cf. p8| . 
As far as the axial perturbation variables are concerned, the set of equations ( 4.1£ ), ( 4.21 ) and ( 4.22| ) makes sense: 
We have three equations governing /iq.;, /ii.z and Ui, for all I. What is not so clear is whether the remaining 0{fl) 
perturbation equations yield unique Ki, Hq.Ii ^2,;, ^Pi and 6ei once the axial variables are specified. (For non-isentropic 
stars, the variables Wi, Vi and Hi i are assumed to be of order il^ so they will not enter into our calculation). In order 
to show that this is the case we must show that the remaining equations can be reduced to five independent ones. 
In this effort we are immediately helped by the fact that, given the assumed ordering (4.18); i) the two equations 
SGtip ~ SnSTt^ and SGte — SirSTtg both imply ( 3.22 ) at order fi, and ii) SGtr — SnSTtr ~ so is automatically 
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satisfied at lower orders. This leaves us with six equations: The equation of state for the perturbations and, for 
example, the five remaining Einstein equations. In other words, the problem would seem to be overdetermined. 
However, for non-isentropic stars the equation of state (2.13) that relates 6p to dp is of order fl^. Thus, we have five 
equations for our five unknown variables and the problem is well-defined. In o ther words, if a discrete mode which 
limits to a purely axial perturbation as 17 — > exists it should follow from ([4.22 ). For completeness, the perturbation 
equations for non-isentropic stars (complete to order il) that follow from ( 4.18 ) are listed in Appendix B. 

Let us now suppose that a distinct r-mode solution exists in the non-isentropic case. One would intuitively expect 
this to be the case since there will then be an internal stratification in the star associated with the entropy gradient. 
In the Newtonian case, this stratification leads to a single r-mode for each combination of I and m at order fl (these 
modes then become non-degenerate at order il^ ||26| , |52| , |53| ) , and it also leads to the presence of non-trivial polar 
g-modes. 



From the above discussion we know that a relativistic r-mode of a non-isentropic star should follow from (4.22). 
We begin our search for such solutions by deriving a constraint on the possible mode-frequencies. We do this by first 
scaling out both f2 and m from the problem by expressing the frequency a in terms of a real parameter a, such that 



Then (4.22) can be written 











dr 





l{l + l) 



1 



2a 



1(1 + 1) 



AM 



hoA > + 16tt{p + p)ahoj = , 



(4.23) 



(4.24) 



where we have used to — Lu/fl. From this equation we see that the eigenvalues a and the corresponding eigenfunctions 
/iQ,; are not explicitly dependent on cither Q or m. The latter means that, if we find an acceptable mode-solution to 



(4.24) it will be relevant for all m ^ for each given multipole I. This would be in accord with the non-isentropic 
Newtonian case where one finds a single r-mode for each combination of I and m at order j2^j5^,|5^ . 

As we will now establish, non-trivial solutions to (4.24) may exist provided that a — Lu vanishes at at least one point 
in the interval r S [0, oo]. To show this we first assume that a — Lj does not have a zero in r G [0, ooj . Then we can 



define a new well-behaved function / through ho^i — r [a — uj)f. By introducing this definition in ( 4.24 ), multiplying 



by f and integrating over r G [0, oo] one can show that (as long as / vanishes both as r and r — > cxj as is 
required by the regularity conditions) 



(a — (D) 



Jo 



{a - Cb)^[l{l + 1) - 2ye^-''\f\^dr . 



(4.25) 



Here both integrands are positive definite, and it follows that we can have no non-trivial solutions for /. 

In other words, a non-trivial solution for /iq,; can only exist if a — a) = at some point in the spacetime. That is, 
the eigenvalue a must lie somewhere in the range 



< a < 1 



(4.26) 



As already noticed by Kojima this agrees well with the Newtonian result. As the star becomes less relativistic 
(I^c — > 1 and our integral relation then predicts a non-trivial solution only for a = 1, i.e. the Newtonian r-mode 
eige nvalue. We will attempt to find discrete r-mode solutions, with frequencies in the permissible range, in Section 

vbI. 



B. The isentropic case 



As indicated above, the conservation of vorticity gives rise to a mixing of axial and polar modes at zeroth order in 
fl. This suggests that the modes of isentropic stars will generically be of a hybrid nature, and as a consequence the 
equations determining the modes are more complicated than those for r-modes of non-isentropic stars. 

The relevant perturbation equations for the isentropic case follow from the spatial components of Eq. ( 4.16| ), which 
for isentropic stars becomes simply, 

Aa;„0 = 0. (4.27) 

We begin by expressing this relation, i.e. 
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= Auj. 



'ap — y a 



A I up 



in terms of ^" and hafs- 



Making use of Eq. (2.11) we have 

A 

where 





e+ p 


Aua 


( n 






n 





A Ur 



(4.28) 



e + p 



Aua = A{gaf3u'^) = Agapuf^ + gapAu" 



(4.29) 



(4.30) 



The ordering ( 4.1^ ) implies that u'^u^hap and g^^hap vanish to zeroth order in fl, since the only zeroth order 
metric components are htr, hte and /ij^. Therefore, 



Au" = u^u^u^Vp^^ 



From Eqs. (2.11) and (3.6) and the relation, 



we obtain, 



to zeroth order in fl. Wc will also use the explicit form of u^p determined from Eq. (fl.2|), 



Tip 



-1/ - 2 • 2 , 

u^p — e ujr sm ( 



and the components of Aua to zeroth order in Q 
Aur — e^^ 

Aue = e'" 



/ 1 \ e2iy 
htr + iK-^ir + f^r^S^ ( —^^ \ + —9^ (r^we-^^^) 

/ite + iK^i^e + VLdei^ - 2w cot 615^ 



Au^ = e 



htip + iK,Q.£^^p + ri9;p^;p + 2[.i; sin 9 cos 6*^9 
+e''9^ {r'^uje-") sin^ 6*6-2^^^ , 



For completeness, we explicitly write down the components of ikJI^ to zeroth order in 17, cf. Eq. ( |4.9| ), 

iKl^e* = 0(^1) = y , ^ , [Vi sin ea^^)™ + mUiY,""] e'"' 

^-^ sm p 



I ^ 



iKn^'p - y , ' . bii-t^rr + t/i sin^a^rr] e^' 

sm 



^-^ sm 



ml g^cri 



inVLS,^ = y i [mViYi" + Ui sinedeY{^] , 



(4.31) 
(4.32) 
(4.33) 
(4.34) 



(4.35) 

(4.36) 
(4.37) 

(4.38) 

(4.39) 
(4.40) 
(4.41) 



(4.42) 
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By making use of Eqs. ( 4.29| ) through (4.41) and the expressions ( 4.42| ) and (O^ ) for the components of inftS^ and 
hai3, we may now write the spatial components of Awa/J- We will use Eq. ( 3.2C ) to eliminate Hi^i (for all I) from the 
resulting expressions and drop the "0" subscript on the metric functions ho.i, writing ho,i = hi. 



dgAu^ ~ du,Aug - dg 



e + p \ e ^ sin ( 



[l{l + l)Kn{hi + Ui)-2mujUi]Y{^ 

I ^ 

-2u;Vi [sinedgY^"' + l{l + 1) cos^F;"] 



^^dr {r^Loe-^") Wi [smOdgY^'^ + 2cos6'Y,"] !>e""* 



(4.43) 



e+p 



dr (e-'Aug) ~ dg (e^'^Au^) 



n J nfl sin 6 



+ ^dr {r^uje-^'') Ui [m2 + l{l + l)(cos2 - 1)] F™ 
-2mdr {uje-^^'Vi) cosOYl^ + ^dr {r^Qe'^") Vi sinSSel^™ 



{AAA) 



Auj^j- — 



e + p 



(e-" Aur) ~ dr {e-" Au^) 



e + p \ e 



Yr 



-2dr {Coe-^^Vi) cos 61 sin 6*99 F;" + {r^ue-^'') ViY{" 

-2mdr {uje-^^'Ui) cos^Y;" + nVldr [e^^^'ihi + C//)] sinedgYC 



(4.45) 



+ ^dr (r^Lue"^'') Ui sin6'(9el^"' U*' 



These equations can be rewritten using the standard identities 

sin^a^rr = IQi+iYi+i - + WiYi-i 
COS 0Yr = Qi+iYr^,+QiYi'I^^ 



(4.46) 
(4.47) 
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with Qi defined as, 



Qi 



{I + m){l — m) 



(2/-l)(2/ + l) 



We then get, from Aive^ = 0, 



= ^1 [1(1 + l)KQ{hi + Ui) - 2mu)Ui] 



'-dr (r^uje^^") Wi + 2{l + l)u;Vi {I - l)QiY{]2i 



Prom Aure = we have, 



+ 



mnndr [e-2-(/i; + Ui)] + 2dr (we-2-[/;) {{I + 1)Q2 _ ;g2^^') 



+^dr (r2we-2-) Ui [w? + l{l + 1) (Q^+i + - l)] 



y-r, 



+ 



QiYi 



1-2 



Prom Aa;,^r = we have. 
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dr [^dr (r2we-2-) Wi] - 2ldr {ue-^^Vi) 



+ 



+ 



Inndr [e-^^ihi + Ui)] - 2mdr {tde-^^Ui) + ^5, (r^we-^-) Ui Qi+iYj™! 

mnndr (e-^^Vi) + 2dr (uje-^^Vi) {{I + l)Qf - IQf^^) 

+ ^dr [r^uie-^-) Vi + dr [Idr [r^Loe-^^] Wi] {Q^^, + - l) 



{l + l)K'[ldr [e-2-(/l; + C/0] 

+2mdr {uje-^^Ui) + ^^9, [r-'ioe-^-) Ui 



+ 



Qi-iQiY{I^2 



dr [^dr (r^we-^-) Wi] + 2{l + l)dr (iJe-^^Vi) 
Finally, let us rewrite these equations using the orthogonality relation for spherical harmonics, 

Y/Yr*dn = Su'Smm', 



where dU is the usual solid angle element on the unit 2-sphere. 
From Auje^p = we have, for all allowed I, 

= [1(1 + l)Kn{hi + Ui) - 2mLdUi] 



+{l + l)Qi ^dr (r^we-^-) Wi-i - 2{l - l)QVi-^ 



~lQi+i 

Similarly, Aw^e = leads to 

0={l-2)Qi_iQi 



^dr 



Wi+i + 2{l + 2)ujVi+i 



+Qi 



-2dr {u>e-^-'Ui-2) + ^^dr (r^we-^-) Ui- 
{l - l)K,ndr (e-^^Vi-i) - 2mdr {uie-^^Vi-i) 



+ 



mnQdr [e-^^ihi + Ui)] + 2dr {uie-^^Ui) {{I + l)Qf - IQf^^) 
+ ^dr (r2a;e-2-) Ui [m^ + l{l + 1) {Q^^, + Of - l) 



-Qi+i 



{I + 2)Kndr {e-^^Vi+i) + 2mdr {luc-^'^Vi+i) 

+=^dr (r^u^e-^'^) y,+r + (/ + 2)K^le-^^ (^^)^ - ^''^1+^ 



+{l + 3)Qi+iQi+2 



2dr {u>e-'-Ui+2) + ^^-^dr {r'ue-''') U1+2 
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and from Aw^r = we get 



+0, 



2{l ~ 2)dr {Qe-^''Vi^2) 



-2mdr (cDe-2-;7,_i) + ^^dr {r^Qe-^-) C/,_i 



(4.55) 



mnUdr {e-'^^'Vi) + 2dr {uoe~'^''Vi) {{I + l)Q} - IQf^^) 

+ ^dr (r^cDe-^-) Vi + dr [Idr {r^Cje~^-) Wi] {Qf^, + Qf - l) 



-Q. 



i+i 



{l + 2)Kndr [e-^^ihi+i + Ui+i)] 



+2mdr {ue-'^'^Ui+i) + 



m{l+2) . 



-Qi+iQi+2 



Wi 



+2 



+ 2{l + 3)dr {uje-^''Vi+2) 



It is instructive to consider the Newtonian limit 



uj{r),iy{r), X{r),hi{r) 0. 



(4.56) 



of these perturbation equations. We have aheady seen that Eq. ( 3.23 ) is the relativistic generahzation of the 
Newt onian mass conservation equation (LF,13) (or Eq. (LF,42)), and that the other 0(1) perturbation equations, 
( |3.20 ) and ( |3.22| ), simply vanish in the Newtonian limit. Similarly, one can readily observe that the conservation of 
vorticity equations have as their Newtonian limit the corresponding equations from Lockitch and Friedman p7| , 



Eq. 
Eq. 
Eq. 



.53) 
,54) 



(|4.55|) 



Eq. (LF,38), 
Eq. (LF,40), 
Eq. (LF,39) 



(and similarly for the other forms of these equations). 

This correspondence leads us to expect the same structure for the relativistic modes as was found in the isentropic 
Newtonian case: we expect to find a discrete set of axial- and polar-led hybrid modes with opposite behavior under 
parity ||2^ . Further, we expect a one-to-one correspondence between these relativistic hybrid modes and the Newtonian 
modes, which the relativistic hybrids should approach in the Newtonian limit. 

In deriving the components of the perturbed vorticity tensor in Newtonian gravity (LF,38)-(LF,40), no assumptions 
about the ordering of the perturbation variables ((5/9, (5w'') in powers of the angular velocity O were required. Thus, 
the theorem concerning the character of the Newtonian modes, [cf. Appendix A, |^] applies to any discrete normal 
mode of a uniformly rotating isentropic star with arbitrary angular velocity. 

We conjecture that the perturbations of relativistic stars obey the same principle: If {£,°',hap) with ^" ^ is 
a discrete normal mode of a uniformly rotating stellar model obeying a one-parameter equation of state, then the 
decomposition of the mode into spherical harmonics Y™ has I = m as the lowest contributing value of /, when m 0; 
and has or 1 as the lowest contributing value of I, when to = 0. 

However, in deriving the curl of the perturbed Eule r equ a tion for relativistic models we have imposed assumptions 
that restrict its generality. We ha ve de rived Eqs. (4.52)-(4.55) in a form that requires a slowly rotating equilib- 
rium model, assumes the ordering ( f.l8 ) and neglects terms of order fl^ and higher. Under these more restrictive 
assumptions, the following theorem holds. 

Theorem 1 Let {ga/3{^),Tap{^)) be a family of stationary, axisymmetric spacetimes describing a sequence of 
stellar models in uniform rotation with angular velocity VL and obeying a one-parameter equation of state, where 
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igaf3(0),Ta/j{0)) is a static spherically symmetric spacetime describing the non-rotating model. Let {^"(fl)^ hapiyt)) 
with ^" ^ &e a family of discrete normal modes of these spacetimes obeying the same one-parameter equation of 
state, where (^"(0), /iq./3(0)) is a stationary non-radial perturbation of the static spherical model. Let (^"(fJo), /ia/3(i^o) 
be a member of this family with D,q <gC ^Ik, the angular velocity of a particle in orbit at the star's equator. Then 
the decomposition o/ (^"(fio), /ia/3(f^o) j'^^o spherical harmonics (i.e., into {l,m) representations of the rotation 
group about its center of mass) has I — m as the lowest contributing value of I, when m ^ 0; and I — 1 as the lowest 
contributing value of I, when m — 0. 

We designate a non-axisymmetric^ mode with parity 
butions only from 



1)™+ an "axial-led hybrid" if and hafs receive contri- 



axial terms with / — m, m + 2, m + 4, ... and 
polar terms with Z = m + 1, m + 3, m + 5, .... 

Similarly, we designate a non-axisymmetricp| mode with parity (—1)™ a "polar-led hybrid" if and h^p receive 
contributions only from 

polar terms with I = m, m + 2, m + 4, ... and 
axial terms with I = m + 1, to + 3, to + 5, .... 

We prove the theorem separately for each parity class in Appendix C. 

In essence, the theorem shows that if a mode of a slowly rotating isentropic star has a stationary non-radial 
perturbation as its spherical limit, then it is generically a hybrid mode with mixed axial and polar angular behavior. 
An immediate consequence of the theorem is that the r-modes of isentropic stars (if they exist at all) must have I — m 
(or Z = 1 if TO = 0) , and it is well known that isentropic Newtonian stars retain a vestigial set of purely axial modes 
- the "classical r-modes" - whose angular behaviour is a purely axial harmonic having I — to. Let us address the 
question of whether or not such pure r-mode solutions also exist in isentropic relativistic stars. 

We apply the perturbation equations for isentropic stars to the case of an axial mode belonging to a pure spherical 
harmonic of index I. In other words, let us assume that hi{r) an d U ijr) (f or so me particular value of I) are the only 
non-vanishing coefficients in the spherical harmonic expansions (4.9) and (4.12) of the Lagrangian displacement, 
and the perturbed metric, ha^, respectively. 

The set of equations to be satisfied are the z er oth o rder (spherical) equations ( 3.2Cl| ), ( 3.22| ) and ( 3.23| ) and the order 
il conservation of circulatio n equ ations ( 1.53 )-( 4.55| ), as well as suitable boundary conditions at infinity and at the 
surface of the star, (Section |IV C ). Recall that as a result of Eq. ( 4.17|), t he co nserva tion of circulation equations are 
linearly dependent and we need only satisfy two of them, say, Eqs. ( 4.53 ) a nd (4.54 ). 

Wi th h i(r) and Ui{r) the only non- vanis hing perturbation variables, Eqs. ( 3.2(]| ) and ( [3.23 ) vanish identically, while 
Eq. (3.22) remains unchanged. Eq. (4.53) becomes, 



= [1(1 + l)Kn{hi + Ui) - 2mLdUi] , 
and Eq. (4.54) with 1^1 + 2,1^1 and 1^1 — 2 gives the equations 



= IQi+iQi^ 



-2dr {Loe-'^''Ui) 



dr {r'^Cje-'^") Ui 



= mnndr [e-^^ihi + Ui)] 

+2dr {Coe-^^Ui) {{I + l)Qf - IQf^^) 

+^dr (r2cSe-2-) Ui [to2 + 1(1 + 1) (Q2^^ +Qf- 1)] , 

I 



(4.57) 

(4.58) 
(4.59) 



= {1 + l)Qi-iQi 



2dr {u 



(4.60) 



®When m = 0, there exists a set of modes with parity -1-1 that may be designated as "axial-led hybrids," since and hafj 
receive contributions only from axial terms with ? = 1, 3, 5, . . . and polar terms with I — 2, 4, 6, . . .. 

^°When m = 0, there exists a set of modes with parity —1 that may be designated as "polar- led hybrids," since and ha/i 
receive contributions only from polar terms with Z = 1, 3, 5, . . . and axial terms with I — 2, 4, 6, . . .. The family of modes 
for which ^" and ha/s receive contributions only from polar terms with Z = 0, 2, 4, . . . and axial terms with Z = 1, 3, 5, . . . would 
have parity -1-1 and could be designated "polar-led hybrids." However, these modes require a more general theorem to establish 
their character. 
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respectively. 

Given the requirement that I = m when m > (and I — 1 when m = 0), one readily finds that these equations can 
be satisfied only when 1 = 1. Thus, no purely axial modes with I = m > 2 exist in isentropic relativistic stars [3^ . The 
dipole {I = 1) solutions turn out to be stationary [a = 0) and have the natural physical interpretation of a uniform 
rotation of the star^. 

In Newtonian isentropic stars there remained a large set of purely axial modes with I = m; the I = m = 2 mode 
being the one expected to dominate the gravitational wave-driven instability of sufficiently hot and rapidly rotating 
neutron stars In isentropic relativistic stars, however, we see that all such pure r-modes with I = m > 2 are 

forbidden by the perturbation equations, and instead must be replaced by a xial-led hybrids. We explicitly construct 



these important hybrid modes to first post-Newtonian order in Section V C 



C. Boundary Conditions 



Having understood the general nature of the relativistic perturbation problem and derived the relevant perturbation 
equations for both isentropic and non-isentropic stars, we want to determine mode-solutions. Before we can do this, 
we need to discuss the boundary conditions that should be imposed. 

For non-isentropic stars, the zeroth order variables are g overned by the sin gle eq uation ( 4.22D , while for isentropic 
stars we have the set of perturbation equations ( 3.2C ), ( 3.22 ), ( 3.23 ) and ( 4.53 )-( |4.55| ). A physically reasonable solution 



(^": hafs) to these equations must be regular everywhere in the spacetime. Of course, the fluid variables Wi{r), Vi{r) 
and Ui{r) (for all I) have support only inside the star, r £ [0, R\. The metric functions Hi i{r ) will also have support 
only inside the star (for all I), since they are directly proportion al to Wi{r) by Eq. (3.20). The metric functions 
hi{r), on the other hand, satisfy a nontrivial differential equation, (3.22), in the exterior spacetime and will, therefore, 
have support on the whole domain r £ [0, oo]. Let us now consider the boundary and matching conditions that our 
solutions must satisfy. 

At the surface of the star, r = R, the perturbed pressure, Ap, must vanish. (Th is is how one defines the surface of 
the perturbed star.) The Lagrangian change in the pressure is given by Eq. (2.11), 



Making use of Eq. ( 1.3(^ ) and the equilibrium equations ( |3.4| ) and (^|^) , we find that at r = i? 



= Ap = 



i?2 (R _ 2Afo 



■^iy/(i?)r,"e*' 



(4.61) 



(4.62) 



where AIq = M{R) is the gravitational mass of the equilibrium star and satisfies 2Mo < R. 

For the equations of state we consideij^ the energy density e(r) either goes to a constant or vanishes at the surface 
of the star in the manner (this would be the behavior for a polytrope). 



e(r) 



(for some constant k). In both cases, it is required that 



R 



Wi{R) = (all I). 



(4.63) 



(4.64) 



If e{R) ^ 0, then Eq. (|I6^ ) requires t his d irectly. On the other han d, if e vanishes as in Eq. (|4.63| ) then Vi{r) will 
diverge at the surface by Eq. ( |3.23 ) if (4.64) is not satisfied. By Eq. (3.20), this also implies that Hi^i{r) vanishes at 
the surface of the star. This boundary condition is (obviously) relevant only in the isentropic case. 



When m = and 1 = 1, the solution corresponds to a small change in the angular velocity of the star about its original 
rotational axis. When I = m = 1, the solution represents uniform rotation of the star about an axis perpendicular to its original 
rotational axis. These solutions are derived in detail by LockitchlSSl and are a generalization of the axial dipole modes studied 
in non-rotating relativistic stars by Campolattaro and Thorne 

^^This restriction can be dropped if the boundary condition Ap{r = i?) = is replaced by Ah{r = R) = 0, with the comoving 
enthalpy h = dp'/{e{p') + p'). 
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In the exterior vacuum spacetime, r > R, we have only to satisfy the single equation ( [3.22| ) for all I, which becomes 



{2-f 



or 



(1 



2Mo 



+ 4Mo 



0, 



(4.65) 



(4.66) 



where we have used e"^"** = (1 — 2Mo/r) for r > R. 

Since this exterior equation does not couple hi{r) for different values of I, we can find its solution explicitly. The 
solution that is regular at spatial infinity can be written 



s=0 



l+s 



(4.67) 



If we substitute this series expansion into Eq. (4.66), we find the following recursion relation for the expansion 
coefficients, 



hi., 



2Mo\ (/ + s-2)(/ + s + l) 



R 



s{2l - 



1) 



(4.68) 



with hi^o an arbitrary normalization constant. We, therefore, have the full solution to zeroth order in of the 
perturbation equations in the exterior spacetime. 

This exterior solution must be matched at the surface of the star to the interior solution for hi{r). One requires 
that the solutions be continuous at the surface. 



lim [hi (i? - e) - hi {R + e)] = 0, 

e— »0 

for all I, and that the Wronskian of the interior and exterior solutions vanish at r = _R, i.e. that 

lim \hi{R ~ e)h'i{R + e) - h[{R - e)hi{R + e)] = 0, 

e— »0 



(4.69) 



(4.70) 



for all I. 

Thus, in solving the perturbation equations to zeroth order in fl we need only work in the interior of the star (as 
i n the Newtonian case). In the interior of a non- isentr opic s tar, the perturbation haf^) must only satisfy Eq. 
( 4.22 ) toge ther wi th th e match ing co nditi o ns (4.69 ) and ( 4.70 ). In the isentropic case we have the full set of c ouple d 
equat ions ( 3.20 ), ( 3.22] ), ( 3.23| ) and ( 4.53| )-(4.55) for all I, subject to the boundary and matching conditions (4.64), 
( p9|) and (FtoI). 

Finally, we note that since we are working in linearized perturbation theory there is a scale invariance to the 
equations. If {S,",ha/3) is a solution to the perturbation equations then [K^" , Khap) is also a solution, for constant 
K. This will not affect the non-isentropic calculation, but in the isentropic case we will find it convenient to impose 
the following normalization condition in addition to the boundary and matching conditions just discussed; 



Um{r = R) = 1 for axial- hybrids, 
Um+i{r = R) = 1 for polar- hybrids. 



(4.71) 



V. RELATIVISTIC CORRECTIONS TO THE R-MODES OF UNIFORM DENSITY STARS 



In a future paper, we will consider the general problem of numerically solving for the r-modes and hybrid modes of 
fully relativistic stars. Preliminary results have already been presented by Lockitch | |3^ . For now, we will focus on 
the post-Newtonian corrections to the Newtonian r-modes. The equilibrium structure of a slowly rotating star with 
uniform density is particularly simple pTf and lends itself readily to such a post-Newtonian analysis. The results we 
obtain in this way provide important insights into the relativistic corrections to the familiar Newtonian r-modes. 
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A. A post-Newtonian uniform density model 



For a spherically symmetric star with constant density, 



e(r) 



3A/o 



(5.1) 



the equilibrium equations (|3.3|)-(3.7) have the well-known exact solution inside the star (r < R), 



p{r) = e 



1 _ 2Mo (7^) 
R \R) 



1 

,21 2 



-2A(r) 



1 



2 1 / _ 2Mo\2 

2 I R 



where A/q is the gravitational mass of the star and R is its radius. 

To find the equilibrium solution corresponding to a slowly rotating star, we must also solve Hartle's 

(PD, 



= r^w" + [4 - r{i^' + X')]ruj' - iriiy' + X')uj 
(see also ^7\) where we may use the spherical solution to write 
r(zy' + A') = 4:TTr^{e+p)e^^ 



(5.2) 

(5.3) 
(5.4) 
(5.5) 

16| equation 
(5.6) 

(5.7) 



2 Ml 



(I_Y 
R \rJ 



2Mo 



.Rj 



To simplify the problem, we expand our equilibrium solution in powers of (2Mq/R) and work only to linear order^. 
We will need the expressions, 



and 



3 1 / r 

2^2 



f2Mo\ , ^ f2Mo 



V R 



V R 



(5.8) 



(5.9) 



Since we are also working to linear order in the star's angular velocity, we may set £7 = 1 without loss of generality. 
We write. 



'i=0 



R 



2i 



(5.10) 



and solve Eq. (5.6) subject to the boundary condition (4.6) at the surface of the star, 



This expansion will give us the first post-Newtonian (IPN) corrections to the Newtonian r- modes. 
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1 = n 



UJ+ -RCu' 



(5.11) 



r=R 



To order (2Mq/R) the solution is, 



""^'^^'-V-^R^ [-IT 



o 



( 2Mo 

V R 



(5.12) 



B. Non-isentropic stars 



In order to find the find the relat ivistic analogue to the familiar Newtonian r-modes of non-isentropic stars, we 
insert the above expressions in (4.22). We also assume that the mode-frequency can be approximated as 



2m 



1(1 + 1) 



and that the eigenfunction takes the form 



1 + Kl 



/2Afo 



V R 



o 



f2Mo 



V R 



(5.13) 



/2Mo 

I R 



o 



/2Mo 

I R 



(5.14) 



Solutions of this form would then lead to Ui ^ 0{1) via (4.21). 

Under these assumptions, (4.22) is trivially satisfied to leading order. At order (2Mo/-R)^ wc find an equation 



1 - 



3r' 
5R^ 



, (0)// _ 1(1 + 



R^ ' 



(5.15) 



Before proceedi ng, it is useful to compare our definition of the post-Newtonian eigenvalue ki to the eigenvalue a we 
used in Section IV A. We then immediately see that 



a = 1 + Ki 



/2Mo 



V R 



and deduce that the established range for possible eigenfrequencics translates into 

- 1 < Ki < . 



(5.16) 



(5.17) 



Within this range there are two possibilities. If ki < —2/5 we will have a singular eigenvalue problem, while for 
—2/5 < Ki < the problem is non-singular. To determine the r-modes to first order in 2Mq/R for the uniform 
density model, one need only consider the simpler non-singular situation, because the eigenvalues of the relativistic 
r-modes turn out to be in the nonsingular range. The continuous part of the spectrum pO| , pl| , as noted earlier, may 
be an artifact of an approximation in which the frequency is real. 
We can rewrite (5.15) as 



1(1 + 1) 



30 



3r2 - 5R^(ki + 1) 



,(0) 







(5.18) 



As long as —2/5 < ki < th is equation can readily be integrated, and the solutions are well-behaved at all values of 
r. We have integrated ( 5.18 ) using a fourth-order Runge-Kutta scheme, initiated from the appropriate regular power 
series solution close to the centre of the star. That is, we use 



Dr 



i+i 



6r' 



R^(Ki + l)[(l + 2)(l + 3) -1(1 + 1)] 



(5.19) 



at an initial point close to r = and then integrate (5.18|) to the surface r = R. At the surface we demand that hf''^ 



and its derivative can be smoothly matched to the exterior solution according to (4.7C). For each value of / we then 
find a single acceptable solution, corresponding to a distinct eigenvalue ki. These eigenvalues, for Z = 2 — 10, are 
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listed in Table |. It should be recalled that the tabulated eigenvalues correspond to mode-frequencies (in the inertial 
frame) given by 



2mn 



l + Kl 



/2Mo 



o 



2Mo 
R 



(5.20) 



A typical eigenfunction, corresponding to Z = 2, is shown in Fig. |^. 

We thus find a single post-Newtonian r-mode solution for each allowed combination of I and m. This is very much 
in accordance with the Newtonian r-mode results for non-isentropic stars at order O (the degeneracy of these modes 
is not broken until at order fi^). The main difference in the relativistic case is that the post-Newtonian corrections 
(of order 2Mq/R) break the degeneracy at order and make it possible for us to determine the eigenf unctions. 

Given these results we expect similar r-mode solutions to exist also in the fully relativistic case. It is, in fact, easy 
to demonstrate this and we have extended our calculation for uniform density star to include all terms in ( 4.22| ) . We 
then find that the m ode-e igenvalue is always such that a — a) = a — cZj/O 7^ in the interior of the star (recall the 
discussion in Section [V A). The solutions to the problem are thus regular. The associated eigenvalues, for Z = 2 — 10 
and a star with compactness M/R = 0.2 are given, and compared to the post-Newtonian results, in Table |. 



C. Isentropic stars 



Having established that discrete r-mode solutions exist for non-isentropic relativistic stars we now turn to the 
isentropic case. As we have shown, we will th en n ot have purel y ax ia l solu tions (for I > 2). Instead, we need to 
calculate hybrid modes by solving Eqs. ( 3.20|) , (3.22), (3.23) and (4.53)-( 4.55| ) subject to the boundary, matching and 
normalization conditions (4.64), ( 4.69| ), ( I.7C| ) and (4.71). As in the non-isentropic case, we seek the post-Newtonian 
corrections to the well-known Newtonian r-modes. For isentropic stars such modes exist only for / = m with frequency 
and radial dependence given by. 



(m+ 1) 

/ 7" \ m+1 

\r) 



(5.21) 
(5.22) 



Therefore, let us make the following Ansatz for o ur pe rturbed solution inside the star. We assume that the coefficients 
of the spherical harmonic expansions (^^) and (4.12) of the Lagrangian displacement, and the perturbed metric, 
ha/3, respectively, have the form. 
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2A/n 



R 
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Wm+l{r) = Wrn.O ( — I (1 



R 
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(5.23) 
(5.24) 
(5.25) 
(5.26) 
(5.27) 
(5.28) 



where ki, Uj^ q, h^.o, hm,i, Wm+i,Oi Vm+i.o, '^m+i.i and Um+2,0 are (as yet) unknown constants. We have chosen 
the form of Um{r) so as to automatically satisfy the normalizati on co ndition ( [4.71 ) and we have chosen the form of 
Wm+i{r) so as to automatically satisfy the boundary condition ( fl.64 ). Note that we have assumed that hi, Vi>, Wu 
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and Ui" are of order (2Mo/i?)^ or higher for all I > m, I' > m + 1 and I" > m + 2. We will justify this Ansatz by 
showing self-consistently that such a s olutio n satisfies the perturbation equations. 

Observe that the exterior solution (4.67) for /im(r) already has a natural expansion in powers of {2AIo/R) as a 
result of the recursion relation ( 4.6g ), 



,RY f2Mo 



O 



2Mo 



(5.29) 



The normalization constant, ft-m.Oj is determined by the matching condition (4.69) 



(5.30) 



while (4.70) imposes the following condition on the interior solution, 



= Kn,o \ -m{h„i,o + hm,i) - [{m + l)h„i,o + {m + 3)Kn,i] 



(5.31) 



or. 



= (2m + l)/i,„.o + {2m + 3)hm,i 



(5.32) 



We turn now to the isentropic perturbation equations, ( 3.20| ), ( |3.22 ), ( 3.23 ) and ( [4.53| )-( 4.55| ). Rec all th at these 
latter three equations are not linearly independent, being related by Eq. ( 4.17] ). Also, because Eq. ( 3.20| ) merely 
expresses Hij{r) in terms oiWi{r)^ we may eliminate Hii{ r) fro m our system and ignore Eq. (3.20). Thus, a complete 
set of perturbation equations is provided by Eqs. (3.22), (3.23), (4.53) an d (J4.54D for all allowed values of /. 



We expand these equations to first pos t-New t onian order using Eqs. (5^), (5.91) and ( 5.12 ) to replace the equilibrium 
quantities and using our Ansatz, Eqs. ( 5.23| )-( 5.2§| ) to replace the various perturbation variables . Th e result is an 
algebraic system of seven independent equations, which, together with our matching condition ( 5.32| ) allows us to 
uniquely find our eight unknown constants, ki, Um,o, hm,o, hm,i, Wm+i,o, Vm+i,o, Vm+i,i and Um+2,o- These equations 
are derived in detail by Lockitch p3|. Here, we will simply present the resulting solution. 



(to + 1) 



4(m- l)(2m+ 11) /2A/o 
5(2m+l)(2m + 5) V~R~ 



O 



/ 2Afo 

V R 



(5.33) 



Urn{r) 



rn+l 



1 + k™.o(1-^) 



o 



/2Mo 

V R 



(5.34) 



(2m + 1) (2m + 3) \R 



2Mo\ , ^ f2Mo 



R 



\ R 



(5.35) 



W^.^ir) ^ (m + l)(m + 2)K^-) ^1 - ^ j ^ ^ j + ^ V i? 



\ /2Afo 



/2Mo 



(5.36) 



/ T- \ m+1 

W(r) = i^(-j 
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where we have defined the constant 



K 



6(m- 1)Q, 



m+l 



5(m + 2)(2m + 5) 



(5.39) 
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and where 



Um,o = - o. / ^^ow'^o' ^ { 48(m + l)4(m + 3)^ (5.40) 



24m(m + 2)(2m + 3) 



+ (2m + 3)2(2m + 5) 



TO(m + 2)2 - 48 



Since our solution satisfies the full perturbation equations to order {2Mq/R), our Ansatz was self-consistent. Thus, 
we have explicitly found the first post-Newtonian corrections to the I = m Newtonian r-modes of isentropic uniform 
density stars. 

The solution reveals the expected mixing of axial and polar terms in the spherical harmonic expansion of . All of 
the isentropic Newtonian r-modes with m > 2 pick up both axial and polar corrections]^ of order (2Mq / R) , becoming 
axial-led hybrid modes of the relativistic star. The I — m = 2 hybrid mode is shown in Figs. ^ and H, and compared 



to the corresponding r-mode in a no n-isentropic star (see Section V B ) . 

In addition, we see from Eq. ( 5.33| ) that the Newtonian r-mode frequency also picks up a small relativistic correction. 



The frequency decreases, just as it does in the non-isentropic case (see Table g), and it is natural that general relativity 
will have such an effect. One reason is that gravitational redshift will tend to decrease the fluid oscillation frequencies 
measured by a distant inertial observer. Also, because these modes are rotationally restored they will be affected by 
the dragging of inertial frames induced by the star's rotation. The Coriolis force is "determined not by the angular 
velocity f2 of the fluid relative to a distant observer but by its angular velocity relative to the local inertial frame, 
Lu{r)." (Hartle and Thorne |5^) Thus, the Coriolis force decreases - and the modes oscillate less rapidly - as the 
dragging of inertial frames becomes more pronounced. 

Finally, we note that the metric perturbation (whose radial behaviour is determined by the function hm) is of 
the same order as the post-Newtonian corrections to the fluid perturbation. Thus, there is no justification for the 
Cowling approximation in constructing the hybrid mode solutions. In Newtonian theory, the Cowling approximation 
corresponds to neglecting the variation in the gravitational potential. The original motivation for this pl[ | is that 
some pulsation modes (in particular the g-modes) are mainly located in the less dense regions close to the surface of 
the star, and do not involve large mass motion. Hence, they will lead to variations in the gravitational potential that 
are small compared to the associated fluid velocities. The obvious generalization of this approach to general relativity 
would be to discard all metric perturbations Jsst . However, as was pointed out by Finn |^6|, this approximation is not 
natural for relativistic g-modes. The main reason is that, even though these modes involve small density perturbations 
they could involve large fluid velocities. Hence, Finn argues that one should keep those metric perturbations that 
can be associated with "momentum transport" in calculations of g-modes. As is easy to see, similar arguments can 
be used for the modes we consider in the present paper. This would suggest that one should not discard the metric 
perturbations hi, Hq and Hi in the relativistic Cowling approximation for r-modes and hybrid modes. Interestingly, 
should we adopt this point of view we retain the main perturbation equations we have used in the present paper. 
Hence, this "approximation" would be consistent with our results. Furthermore, this would explain why the attempts 
to find relativistic r-modes within the Cowling approximation (by neglecting all metric perturbations) have failed [ p2| . 
Of course, this discussion has little relevance for the present study. But it could be of crucial importance for attempts 
to find r-modes in numerical simulations (by studying fluid motion in relativistic simulations with a "frozen" metric) 
that are currently under way |^,^^. 

VI. DISCUSSION 

In this paper we have taken the first steps towards an understanding of both r-modcs and rotational hybrid modes 
of rotating relativistic stars. We have derived the perturbation equations that govern these modes to linear order in 
the rotation frequency f2 (at which the star is still spherical). For non-isentropic stars we have focused on modes that 
have a purely axial limit as 17 — > 0. These would be a natural relativistic generalization of the Newtonian r-modes. 
For isentropic stars (and multipoles / > 2) we have shown that no such modes exist in the relativistic case, even 
though Newtonian stars retain a vestigial set corresponding to Z = m. Instead, all modes of isentropic stars must have 
a hybrid nature. Having derived the relevant perturbation equations we calculate relativistic corrections at the first 
post-Newtonian level (order 2Mo/R) to the Newtonian r-modes of both non-isentropic and isentropic stars. 



When m = 1, the constant K vanishes and we recover the axial dipole solution mentioned in Section IV B 
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It is worth pointing out that, even though our results for isentropic and non-isentropic stars are quite different, the 
particular modes that we have focussed on (the analogues of the vestigial I = m r-modes that remain for isentropic 
Newtonian stars) are not too dissimilar. As is clear from the results given in Table |, the mode frequencies in the 
two cases we have considered do not differ by more than a few percent. Furthermore, we can see from Figure ^ that 
the axial eigenfunctions /iq,; are similar. There are, of course, still considerable differences between the two cases. In 
the non-iscntropic case we predict that purely axial modes exist for all combinations of I and m ^ 0, while in the 
isentropic case all modes are hybrids. Still, the fact that our results for the two cases seem consistent is encouraging. 
We anticipate that further work will eventually unveil a behavior quite similar to that of the Newtonian problem, for 
which the detailed isentropic limit has been investigated by Yoshida and Lee . 

This paper represents progress in several important directions, but a considerable amount of work remains before 
we can claim to have a complete understanding of the nature of the r-modes and hybrid modes in relativity. For 
example, we have not yet discussed how the inferred changes in both mode-frequency and eigenfunction will affect the 
strength of the gravitational-wave driven instability. To do this we need to estimate the rate at which these modes 
radiate gravitational waves, and also assess the strength of various dissipation mechanisms (like viscosity) that tend 
to damp an unstable mode. This is obviously an important issue and we plan to address it once our ongoing work 
on fully relativistic hybrid modes of isentropic stars is completed. At that point it will also be appropriate to obtain 
and discuss results for different realistic equations of state. 
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APPENDIX A: EQUATIONS DESCRIBING STATIONARY PERTURBATIONS OF SPHERICAL STARS 

We have derived the various equatio ns go v ernin g stationary per turbat ions o f a spherical star using the Maple tensor 
package by substituting expressions ( ^.ll| )-( ^?T^ ) into Eqs. (3^) and ( 3.10|) (making hberal use of the equihbrium 
equations (3.4) through (3.7) to simplify the expressions). The resulting equations are hsted in the three distinct 
cases I > 2, I = 1 and / = below. 



1. The case I > 2 



The non- vanishing components of the perturbed Einstein equation for I > 2 are as follows. We will use Eq. (A4) 
below, to replace H2 by Hq. From SG/ ~ SnSTf.* we have (using primes to denote derivatives with respect to r) 



= e-^\^K" + e-2^(3 - rX')rK' - + 1) - l] K 



(Al) 



From 5G^^ ~ SnST,,^ we similarly have 

= e-2^(l + rv')rK' - + l)-l\K 



-2A 



riJ^ + + 1) - 1 - Sttt^p] Ho - Sirr^Sp. 
From + 5G^ = Stt ((5T/ -t- we have 

= e~2A^2^// _^ g-2A ^^(^^1 _ _^ 2] rK' - IQirr'^Sp 

-e-^^r^iJo - e"2^(3riy' - r\' + 2)rH'^ - mnr^pHo- 
From SG/ - 5G^ = Stt ((ST/ - 5T/) we have 



(A2) 



(A3) 
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H2 — Hq. 



From SGf = SttSTJ we have 



From 5G/ = 'SiTtSTJ we have 



From = Stt^Tj ^ we have 



From JGj'^ = ^^T5T^^ we have 



+ 1) 



= e-^"'-^) 



e("'-A)ijJ + 167r(e + p)e2V. 



From = STr^r^"^ we have 
Finally, from 5GJf = ^TiSTg^ we have 



-e {v + \ ) 2 



/io = -(i^' + A')C/. 
r 



(l^l){l + 2)hi =0. 



(A4) 
(A5) 

(A6) 

(A7) 

(A8) 
(A9) 
(AlO) 



= e^^^rH^ + (2 - 87rr^e) iJa - STrr^Se. 



From (5G/ = Stt^TL'' we have 



2. The case / = 1 

The / = 1 case differs from / > 2 in two respects Firstly, H2{r) ^ Ho{r), because the equation 5Gg^ — SG^ = 
8tt {5Tg^ — 5T^^ vanishes identically. Secondly, we may exploit the aforementioned gauge freedom for this case to 
eliminate the metric functions K{r) and /ii(r). (We note that Eq. (A£) implies hi{r) = for I > 2 anyway.) With 
these two differences taken into account the non-vanishing components of the perturbed Einstein equation for / = 1 
are as follows. From (5Gj* = SttST^ * we have 

(All) 
(A12) 

(A13) 
(A14) 
(A15) 

(A16) 
(A17) 



= e^^^rHo ~Ho + {l + &T:r^p) H2 + Snr^Sp. 
From 5Gf + 6G/ = 87r ((5T/ + (JT/) we have 

= e^2A^2^^/ ^ e-^^{2rv' - rX' + l)rHo - Hq 

+e"^^(l + ri^')rH^ + (1 + 16nr'^p)H2 + IGnr^Sp 



From SG/ = Stt^T/ we have 
From SG/ — SttJTj'' we again have 
From (5G/ = SttST^ ^ we again have 



Q = rH'a + {rv' - l)Ho + [rv' + 1)^2 
Hi + 8'K{e+p)e^^rW. 



+ l&TT{€+p)e^^V. 



Finally, from 5Gi^ = 8^T5T^ we have 



2 2 



/lo = -(j^' + A')C/. 
r 
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3. The case / = 



The I — case differs yet again from the previous two, being the case of stationary, sphericaUy symmetric per- 
turbations of a static, spherical equihbrium. To maximize the similarity to the preceding two cases we will use the 
same form for the perturbed metric except that we may now exploit the gauge freedom for this case to eliminate the 
functions K{r), Hi(r) and hi{r). The non- vanishing components of the perturbed Einstein equation for I = are as 
follows. 

From ^G/ = STr^r^ * we have 

= e^^^rH^ + (l - 87rr^e) H2 - Snr^Se. (A18) 

From = 8nST/ we have 

= e-^^rH'o + (l + Sirr^p) H2 + Snr^Sp. (A19) 
From 6G/ + = Stt {5T/ + 5T/) we have 

= e-^^r^Ho + e-^^{2rv' - rX' + l)rHo 

+e"^^(l + ri^')rH2 + l&T^r^pH^ + Unr^dp (A20) 

Finally, from SGf^ = SttST^^ we have 

= 167r(e + p)W. (A21) 



APPENDIX B: PERTURBATION EQUATIONS FOR SLOWLY ROTATING NON-ISENTROPIC STARS 



The assumption of a purely axial perturbation in the limit 51 leads to the following equations (cf. ( 4.18| )): The 
three axial quantities follow from 



(7 + mf2 



2mLU 





xd_ 




\dho,i 






dr 
= . 


dr 




+ mfl 


2m(jj 


Ui + { 


(7 + r 


l{l 


+ 1)J 



1(1 + 1) 4M 



ho,i 



and 



l{l + 1) > i{o' + muj)e 



ho., 



{l-l)il + 2)hi^i 



— 2imuj'e ^"^hoj = , 



(Bl) 
(B2) 

(B3) 



or, alternatively, equation ( 4.2Cl| ). 

The solutions to these three equations then serve as sources for three of the remaining (at order fl) Einstein 
equations that determine the polar parity metric perturbations: 

{I - + 1){1 + 2)e'^H2j - Ho,i)Yr 

- {r^e-^^Lo'h'o^i + [l{l + l)uj - 2re-^^Lu' - Wirr^ip + e)Lu] hoj - 16TTi{p + eyCjUi} 

X {2{l + 2) smedeYi"" + + 1) cos6'r;'"} = (B4) 



1(1 + 1)6^" [r{K'i - H'^^i) + (1 - rv') Ho, - (1 + ri,') H^j] 1^" 

+ 1) {2rujh'Q i + [ruj' - 2lu{1 + rv')]ho,i} shiOdeYi"^ 
+2ruj'hQj{smedgYi"' - l{l + 1) cos^Y,™) = 



(B5) 



and 
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1(1 + 1)6 



2v 



-{Ki-Ho,i) + 
il-l)il + 2)r + 6M 



6M 



- Sire HoA 



sinedoY,''' 



^uj'e-^^ho,i - S2TTi{p + e)QUi ) [sin edgYi"" - 1(1 + 1) cos OYi"') 



Luhoicoser,"' = 



(B6) 



These equations determine the polar metric perturbations Ki, Hq^i and i?2,i once /iq,; and Ui are known for aU I. 
FinaUy, the last two Einstein equations lead to the following equations for Spi and Sei (recall that for non-iscntropic 
stars the equation of state does not link these two quantities until at order fl^); 



— + 2uj'- 2uov'] h'a I + {^Loiv'f - 2oj'v' - —v' ] ho,i + 2,2m{p + ejQe^^Ui 



smOdgY/' 



-^fc^c^e^^V; cos^y/" = 



(B7) 



AiT{p + e)re'^\K[ + -{Hii-Hii) 



^ ^ -e2^(i7o,i + H2.1) + Snip + e)e^^H2,i + Sne'^^Spi + Sei) 



2r2 



Yr 



—oje 
r 



-2uut 



h'^ i + (^{uj' - 2ujv') + lQiTe^^{p + e)u} - 
+167ri(p + €)oje^^~^''Ui] sin 6*96 F/" 



'ho.icoseY,"' = 



(B8) 



APPENDIX C: PROOF OF THEOREM |l| 
1. Axial-led hybrids with m > 0. 



Let I be the smallest value of I' for which Ui' ^ in the spherical harmon ic exp ansion (O) of the displacement 
vector or for which hii = ho^i' 7^ in the spherical harmonic expansion ( [4.12[ ) of the metric perturbation hafj. 
The axial parity of (^", hafj), (—1)'+^, and the vanishing of YJ'" for I < m implies / > m. That the mode is axial-led 
means Wi' = 0, Vii = and Hi 11 — f or I' < I. We show by contradiction that I = m. 

Suppose l>m+l. From Eq. ( [4.53| ), / AuJevY{"'^dft = 0, we have 



1(1 + l)Kn{hi + Ui) - 2mQUi 



(CI) 



-IQ 



-dr {r^uie-^") Wi+i + 2{l + 2)ujVi+i 
and from Eq. (4.55) with I ^ I — I, J Auj^rYi*^d^ — 0, we have 

= -{{1 + l)Kndr [e-^^ihi + Ui)] + 2mdr (ioe-^^Ui) + ^^^^9^ pcSe-^"') U, 



(C2) 



'Qi+i 



dr [Idr (r'^uje' 



2v\ 



Wi 



+1 



2{l + 2)dr {Coe-^^Vi 



+1; 
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Together these give, 



2 {r^ue-^'') ' dr 



or, 



Ui = K {oje'^'') 2 
(for some constant K) which is singular as r — > 0. 



('+2) 

r 



(C3) 
(C4) 

(C5) 



2. Axial-led hybrids with m = 0. 



Let m = and let I be the smallest value of I' for which Ui' 7^ in the spherical harmoni c exp ansion (4.9) of 
the displacement vector or for which /i;/ = ho.i' 7^ in the spherical harmonic expansion ( 4.12 ) of the metric 
perturbation hap. Since VaYo — 0, the mode vanishes unless I > 1. That the mode is axial-led means Wi> = 0, 
Vi' = and -ffi.c = for I' < I. W e show by contradiction that I — 1. 

Suppose I > 2. Then Eq. ( [4.54D with I ^ I - 2, J At^^ei^l^aC^f^ = 0, becomes 



= 2dr {Coe-^^Ui) + —dr (r^we-^") Ui 



2 (r^bje 



r \ijje ) ^ Ui 



or, 



Ui = K {oje-^'') 2 
(for some constant K) which is singular as r — ^ 0. 



('+2) 

r 



(C6) 
(C7) 

(C8) 



3. Polar-led hybrids with m > 0. 



Let I be the smallest value of /' for which Wf 7^ or V/' 7^ in the spherical harmonic expansion (49) of the 
displacement vector f", or for which Hi i' =/= in the spherical harmonic expansion ( 4.12| ) of the metric perturbation 
hap. The polar parity of {^^^jha/i), (—1)', and the vanishing of YJ'" for Z < m implies I > m. That the mode is 
polar-led means Ui' = and hii = for /' < /. Wc show by contradiction that I — m when m > and that I — 1 
when m — 0. 



Suppose / > TO + 1. From Eq. ( [l.53[ ) with / / - 1, / /^Ldg^Yi^dVl = 0, we have 



= (/ - l)Qi 



dr {r^uie-^") Wi + 2{l + l)QVi 



(C9) 



Substituting for Vi using Eq. (B.23), we find 



= -dr (r^Cde-^") Wi + 2oje-^''- -dr 

r ' (e+P) 



2 r^oje 



2-,„-2i/\-2"~ Q 



(CIO) 
(Cll) 



with solution, 
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Wi=K iCoe-^")"^ -r-C+i) (C12) 

(for some constant K) which is singular as r ^ 0. 

When ni — this argument fails to establish that / cannot be equal to 1, because Eq. ( |C^ ) is trivially satisfied for 
Z = 1 as a result of the overall / — 1 factor. Instead, the argument proves that I cannot be greater than 1 in this case 
and therefore that 1 = 1. 
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0.04 




FIG. 1. Numerically determined post-Newtonian r-modo cigenfunction hi{r) for Z = 2 in the interior of a uniform density star 
(solid line). The result is compared to the corresponding cigenfunction for the particular hybrid mode that is the relativistic 
counterpart of the Newtonian I = rn — 2 r-mode in an isentropic star (shown as a dashed curve). Of course, in the isentropic 
case several other functions (such as Wm+i and Vm+i) are also non-zero (see Fig. 2). The functions are normalized to agree at 
the surface of the star and the overall scale is set by the normalization of Fig. 2. 
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FIG. 2. The (r/R)^ radial dependence of the Newtonian I — m = 2 r-mode is shown (dashed curve) together with the 
post-Newtonian corrections to this mode for a uniform densit y sta r of compactness 2M/R = 0.2, i.e., the coefficients Ui(r), 
Wiir), and Vi{r) with Z < 4 of the spherical harmonic expansion (4.9) (solid curves). The vertical scale is set by the normalization 
of U2(r) to unity at the surface of the star, and the other coefficients have been scaled by a factor of 100. Thus, while the 
relativistic corrections to the equilibrium structure of the star are of order 20%, the relativistic corrections to the r-mode are 
only of order 1%. 
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TABLE I. Relativistic r-mode and hybrid mode frequencies of uniform density stars. We list the numerically determined 
values of the post-Newtonian r-mode frequency correction, ti^ ^^^^ , for the non-isentropic star and compare the corresponding 

eigenvalues Ofpjv to ones obtained for fully relativistic uniform density stars, a, for a star of compactness 2Mq/R = 0.4. In this 
case the value of the frame dragging at the surface of the star leads to Lj{R)/i} = 0.84424 and we can see that the eigenvalues 
approach this value as I increases. It is also interesting to compare our post-Newtonian eigenvalues to the result we deduce 
for the hybrid I — m modes of isentropic stars, '^ij^yj^j.jjj = — 4(m — l)(2m + ll)/5(2m + l)(2m + 5). The two results typically 
do not differ by more than a few percent. This is important since the two modes (for / = rn) correspond to the relativistic 
analogue (for uou-isciitropic and isentropic stars, rt'spccth'cly) of tfic saiu(> .Xtnvtouiaii r-modc. 



/ 


'^^hybrid 


r-modc 


CtpN 


a 


2 


-0.2667 


-0.2629 


0.8949 


0.9086 


3 


-0.3532 


-0.3428 


0.8629 


0.8699 


4 


-0.3897 


-0.3734 


0.8506 


0.8561 


5 


-0.4073 


-0.3868 


0.8453 


0.8502 


6 


-0.4163 


-0.3931 


0.8428 


0.8474 


7 


-0.4211 


-0.3962 


0.8415 


0.8460 


8 


-0.4235 


-0.3979 


0.8408 


0.8453 


9 


-0.4247 


-0.3988 


0.8405 


0.8448 


10 


-0.4251 


-0.3993 


0.8403 


0.8446 
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